Transcriptomic changes in the microsporidia proliferation and host responses in congenitally infected embryos and larvae.

Congenital infection caused by vertical transmission of microsporidia N. bombycis can result in severe economic losses in the silkworm-rearing industry. Whole-transcriptome analyses have revealed non-coding RNAs and their regulatory networks in N. bombycis infected embryos and larvae. However, transcriptomic changes in the microsporidia proliferation and host responses in congenitally infected embryos and larvae remains unclear. Here, we simultaneously compared the transcriptomes of N. bombycis and its host B. mori embryos of 5-day and larvae of 1-, 5- and 10-day during congenital infection. For the transcriptome of N. bombycis, a comparison of parasite expression patterns between congenital-infected embryos and larva showed most genes related to parasite central carbon metabolism were down-regulated in larvae during infection, whereas the majority of genes involved in parasite proliferation and growth were up-regulated. Interestingly, a large number of distinct or shared differentially expressed genes (DEGs) were revealed by the Venn diagram and heat map, many of them were connected to infection related factors such as Ricin B lectin, spore wall protein, polar tube protein, and polysaccharide deacetylase. For the transcriptome of B. mori infected with N. bombycis, beyond numerous DEGs related to DNA replication and repair, mRNA surveillance pathway, RNA transport, protein biosynthesis, and proteolysis, with the progression of infection, a large number of DEGs related to immune and infection pathways, including phagocytosis, apoptosis, TNF, Toll-like receptor, NF-kappa B, Fc epsilon RI, and some diseases, were successively identified. In contrast, most genes associated with the insulin signaling pathway, 2-oxacarboxylic acid metabolism, amino acid biosynthesis, and lipid metabolisms were up-regulated in larvae compared to those in embryos. Furthermore, dozens of distinct and three shared DEGs that were involved in the epigenetic regulations, such as polycomb, histone-lysine-specific demethylases, and histone-lysine-N-methyltransferases, were identified via the Venn diagram and heat maps. Notably, many DEGs of host and parasite associated with lipid-related metabolisms were verified by RT-qPCR. Taken together, simultaneous transcriptomic analyses of both host and parasite genes lead to a better understanding of changes in the microsporidia proliferation and host responses in embryos and larvae in N. bombycis congenital infection.

The BTB-ZF gene Bm-mamo regulates pigmentation in silkworm caterpillars.

The color pattern of insects is one of the most diverse adaptive evolutionary phenotypes. However, the molecular regulation of this color pattern is not fully understood. In this study, we found that the transcription factor Bm-mamo is responsible for black dilute (bd) allele mutations in the silkworm. Bm-mamo belongs to the BTB zinc finger family and is orthologous to mamo in Drosophila melanogaster. This gene has a conserved function in gamete production in Drosophila and silkworms and has evolved a pleiotropic function in the regulation of color patterns in caterpillars. Using RNAi and clustered regularly interspaced short palindromic repeats (CRISPR) technology, we showed that Bm-mamo is a repressor of dark melanin patterns in the larval epidermis. Using in vitro binding assays and gene expression profiling in wild-type and mutant larvae, we also showed that Bm-mamo likely regulates the expression of related pigment synthesis and cuticular protein genes in a coordinated manner to mediate its role in color pattern formation. This mechanism is consistent with the dual role of this transcription factor in regulating both the structure and shape of the cuticle and the pigments that are embedded within it. This study provides new insight into the regulation of color patterns as well as into the construction of more complex epidermal features in some insects.

Meiosis: The silk moth and the elephant.

In most eukaryotes, balanced chromosome segregation at meiosis requires crossovers, but female Bombyx mori lack these structures. Instead, the synaptonemal complex is repurposed to compensate for this absence of crossovers, a remarkable example of exaptation.

Lipid homeostasis is essential for oogenesis and embryogenesis in the silkworm, Bombyx mori.

Reproduction, a fundamental feature of all known life, closely correlates with energy homeostasis. The control of synthesizing and mobilizing lipids are dynamic and well-organized processes to distribute lipid resources across tissues or generations. However, how lipid homeostasis is precisely coordinated during insect reproductive development is poorly understood. Here we describe the relations between energy metabolism and reproduction in the silkworm, Bombyx mori, a lepidopteran model insect, by using CRISPR/Cas9-mediated mutation analysis and comprehensively functional investigation on two major lipid lipases of Brummer (BmBmm) and hormone-sensitive lipase (BmHsl), and the sterol regulatory element binding protein (BmSrebp). BmBmm is a crucial regulator of lipolysis to maintain female fecundity by regulating the triglyceride (TG) storage among the midgut, the fat body, and the ovary. Lipidomics analysis reveals that defective lipolysis of females influences the composition of TG and other membrane lipids in the BmBmm mutant embryos. In contrast, BmHsl mediates embryonic development by controlling sterol metabolism rather than TG metabolism. Transcriptome analysis unveils that BmBmm deficiency significantly improves the expression of lipid synthesis-related genes including BmSrebp in the fat body. Subsequently, we identify BmSrebp as a key regulator of lipid accumulation in oocytes, which promotes oogenesis and cooperates with BmBmm to support the metabolic requirements of oocyte production. In summary, lipid homeostasis plays a vital role in supporting female reproductive success in silkworms.

Harvesting Silk Fibers from Plodia interpunctella: Role of Environmental Rearing Conditions in Fiber Production and Properties.

Silk fibers are produced by a wide variety of insects. The silkworm Bombyx mori (Bombyx) was domesticated because the physical properties of its silk fibers were amenable to the production of fine textiles. Subsequently, engineers have regenerated silk fibroin to form biomaterials. The monocular focus on Bombyx silk has underutilized the expanse of diverse silk proteins produced by more than 100,000 other arthropods. This vast array of silk fibers could be utilized for biomedical engineering challenges if sufficient rearing and purification processes are developed. Herein, we show that the moth, Plodia interpunctella (Plodia), represents an alternative silk source that is easily reared in highly regulated culture environments allowing for greater consistency in the silk produced. We controlled the temperature, resource availability (larvae/gram diet), and population density (larvae/mL) with the goal of increasing silk fiber production and improving homogeneity in Plodia silk proteins. We determined that higher temperatures accelerated insect growth and reduced life cycle length. Furthermore, we established initial protocols for the production of Plodia silk with optimal silk production occurring at 24 °C, with a resource availability of 10 larvae/gram and a population density of 0.72 larvae/mL. Population density was shown to be the most prominent driving force of Plodia silk mat formation among the three parameters assessed. Future work will need to link gene expression, protein production and purification, and resulting mechanical properties as a function of environmental cues to further transition Plodia silk into regenerated silk fibroin biomaterials.

Effective expression and characterization of the receptor binding domains in SARS-CoV-2 Spike proteins from original strain and variants of concern using Bombyx mori nucleopolyhedrovirus in silkworm.

A new coronavirus, known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is responsible for the global pandemic of COVID-19 in 2020. Through structural analysis, it was found that several amino acid residues in the human angiotensin-converting enzyme-2 (hACE2) receptor directly interact with those in the receptor binding domain (RBD) of the spike glycoprotein (S-protein). Various cell lines, including HEK293, HeLa cells, and the baculovirus expression vector system (BEVS) with the insect cell line Sf9, have been utilized to produce the RBD. In this study, we investigated the use of Bombyx mori nucleopolyhedrovirus (BmNPV) and BEVS. For efficient production of a highly pure recombinant RBD protein, we designed it with two tags (His tag and STREP tag) at the C-terminus and a solubilizing tag (SUMO) at the N-terminus. After expressing the protein using BmNPV and silkworm and purifying it with a HisTrap excel column, the eluted protein was digested with SUMO protease and further purified using a Strep-Tactin Superflow column. As a result, we obtained the RBD as a monomer with a yield of 2.6 mg/10 mL serum (equivalent to 30 silkworms). The RBD showed an affinity for the hACE2 receptor. Additionally, the RBDs from the Alpha, Beta, Gamma, Delta, and Omicron variants were expressed and purified using the same protocol. It was found that the RBD from the Alpha, Beta, Gamma, and Delta variants could be obtained with yields of 1.4-2.6 mg/10 mL serum and had an affinity to the hACE2 receptor.

Comparative analysis of LC-MS/MS and real-time PCR assays for efficient detection of potential allergenic silkworm.

The silkworm (Bombyx mori) has long been valued food and feed in East Asia for its abundant nutritional and medicinal attributes, conversely, it can elicit allergic responses in susceptible individuals. Therefore, the development of silkworm detection method is required to avert allergenic incidents. In this study, two methodologies, tandem mass spectrometry (LC-MS/MS) and real-time PCR, were developed to achieve effective silkworm detection. These methods exhibited exceptional sensitivity in identifying silkworm presence in processed foods. Furthermore, model cookies spiked with silkworm were used to validate the sensitivities of LC-MS/MS (0.0005%) and real-time PCR (0.001%). Overall, these techniques were useful for trace silkworm detection in food products; therefore, they may help prevent allergic reactions. To the best of our knowledge, this study represents the first comparison of LC-MS/MS and real-time PCR methods for silkworm detection, marking an important contribution to the field. Data are available from ProteomeXchange under identifier PXD042494.

Bmpali, Bmb1 and Bmcap are necessary for uric acid granule formation in Bombyx mori.

Uric acid is the end-product of nitrogen metabolism of the silkworm and other lepidopterans. The accumulation of uric acid particles in the epidermis causes the larval silkworm to appear white and opaque. However, the mechanism of uric acid granule formation is still unclear. Silkworm epidermis color is linked to the genes responsible for uric acid particle formation. We first identified two genes in the Bombyx mori genome that encode subunits of the Bloc-1 (Biogenesis of Lysosome-related Organelles Complex-1) by homology to these genes in other eukaryotes, Bmpali and Bmb1. Mutation in these genes caused a transparent phenotype in the silkworm larvae, and the loss of BmBloc-1 subunit gene Bmcap resulted in the same phenotype. These three genes are highly conserved between human and silkworm. We discovered that Bmpali, Bmcap, and Bmb1 localize in the cytoplasm of BmN cells. Yeast two-hybrid assays demonstrated that the Bmpali physically interacts with both Bmcap and Bmb1. Investigating the roles of Bmpali, Bmb1, and Bmcap is essential for uric acid granule formation understanding in Bombyx mori. These mutants present a valuable silkworm model for studying the biogenesis of lysosome-related organelles (LROs).

Structural Insights into the Substrate Binding of Farnesyl Diphosphate Synthase FPPS1 from Silkworm, Bombyx mori.

Farnesyl diphosphate synthase (FPPS) is an important enzyme involved in the juvenile hormone (JH) biosynthesis pathway. Herein, we report the crystal structure of a type-I Lepidopteran FPPS from Bombyx mori (BmFPPS1) at 2.80 Å resolution. BmFPPS1 adopts an α-helix structure with a deep cavity at the center of the overall structure. Computational simulations combined with biochemical analysis allowed us to define the binding mode of BmFPPS1 to its substrates. Structural comparison revealed that BmFPPS1 adopts a structural pattern similar to that of type-II FPPS but exhibits a distinct substrate-binding site. These findings provide a structural basis for understanding substrate preferences and designing FPPS inhibitors. Furthermore, the expression profiles and RNA interference of BmFPPSs indicated that they play critical roles in the JH biosynthesis and larval-pupal metamorphosis. These findings enhance our understanding of the structural features of type-I Lepidopteran FPPS while providing direct evidence for the physiological role of BmFPPSs in silkworm development.

A major endogenous glycoside hydrolase mediating quercetin uptake in Bombyx mori.

Quercetin is a common plant flavonoid which is involved in herbivore-plant interactions. Mulberry silkworms (domestic silkworm, Bombyx mori, and wild silkworm, Bombyx mandarina) take up quercetin from mulberry leaves and accumulate the metabolites in the cocoon, thereby improving its protective properties. Here we identified a glycoside hydrolase, named glycoside hydrolase family 1 group G 5 (GH1G5), which is expressed in the midgut and is involved in quercetin metabolism in the domestic silkworm. Our results suggest that this enzyme mediates quercetin uptake by deglycosylating the three primary quercetin glycosides present in mulberry leaf: rutin, quercetin-3-O-malonylglucoside, and quercetin-3-O-glucoside. Despite being located in an unstable genomic region that has undergone frequent structural changes in the evolution of Lepidoptera, GH1G5 has retained its hydrolytic activity, suggesting quercetin uptake has adaptive significance for mulberry silkworms. GH1G5 is also important in breeding: defective mutations which result in discoloration of the cocoon and increased silk yield are homozygously conserved in 27 of the 32 Japanese white-cocoon domestic silkworm strains and 12 of the 30 Chinese ones we investigated.

K64 acetylation of heat shock protein 90 suppresses nucleopolyhedrovirus replication in Bombyx mori.

HSP90 is a highly conserved chaperone that facilitates the proliferation of many viruses, including silkworm (bombyx mori) nucleopolyhedrovirus (BmNPV), but the underlying regulatory mechanism was unclear. We found that suppression of HSP90 by 17-AAG, a HSP90-specific inhibitor, significantly reduced the expression of BmNPV capsid protein gp64 and viral genome replication, whereas overexpression of B. mori HSP90(BmHSP90) promoted BmNPV replication. Furthermore, in a recent study of the lysine acetylome of B. mori infected with BmNPV, we focused on the reduced viral proliferation due to changes of BmHSP90 lysine acetylation. Site-directed introduction of acetylated (K/Q) or deacetylated (K/R) mimic mutations into BmHSP90 revealed that lysine 64 (K64) acetylation activated the JAK/STAT pathway and reduced BmHSP90 ATPase activity, leading to diminished chaperone activity and ultimately inhibiting BmNPV proliferation. In this study, a single lysine 64 acetylation change of BmHSP90 was elucidated as a model of posttranslational modifications occurring in the wake of host-virus interactions, providing novel insights into potential antiviral strategies.

Exploring the apoptotic effects of sericin on HCT116 cells through comprehensive nanostring transcriptomics and proteomics analysis.

Sericin, a silk protein from Bombyx mori (silkworms), has many applications, including cosmetics, anti-inflammation, and anti-cancer. Sericin complexes with nanoparticles have shown promise for breast cancer cell lines. Apoptosis, a programmed cell death mechanism, stops cancer cell growth. This study found that Sericin urea extract significantly affected HCT116 cell viability (IC50 = 42.00 ± 0.002 µg/mL) and caused apoptosis in over 80% of treated cells. S-FTIR analysis showed significant changes in Sericin-treated cells' macromolecule composition, particularly in the lipid and nucleic acid areas, indicating major cellular modifications. A transcriptomics study found upregulation of the apoptotic signaling genes FASLG, TNFSF10, CASP3, CASP7, CASP8, and CASP10. Early apoptotic proteins also showed that BAD, AKT, CASP9, p53, and CASP8 were significantly upregulated. A proteomics study illuminated Sericin-treated cells' altered protein patterns. Our results show that Sericin activated the extrinsic apoptosis pathway via the caspase cascade (CASP8/10 and CASP3/7) and the death receptor pathway, involving TNFSF10 or FASLG, in HCT116 cells. Upregulation of p53 increases CASP8, which activates CASP3 and causes HCT116 cell death. This multi-omics study illuminates the molecular mechanisms of Sericin-induced apoptosis, sheds light on its potential cancer treatment applications, and helps us understand the complex relationship between silk-derived proteins and cellular processes.

High-throughput and genome-scale targeted mutagenesis using CRISPR in a nonmodel multicellular organism, Bombyx mori.

Large-scale genetic mutant libraries are powerful approaches to interrogating genotype-phenotype correlations and identifying genes responsible for certain environmental stimuli, both of which are the central goal of life science study. We produced the first large-scale CRISPR-Cas9-induced library in a nonmodel multicellular organism, Bombyx mori We developed a piggyBac-delivered binary genome editing strategy, which can simultaneously meet the requirements of mixed microinjection, efficient multipurpose genetic operation, and preservation of growth-defect lines. We constructed a single-guide RNA (sgRNA) plasmid library containing 92,917 sgRNAs targeting promoters and exons of 14,645 protein-coding genes, established 1726 transgenic sgRNA lines following microinjection of 66,650 embryos, and generated 300 mutant lines with diverse phenotypic changes. Phenomic characterization of mutant lines identified a large set of genes responsible for visual phenotypic or economically valuable trait changes. Next, we performed pooled context-specific positive screens for tolerance to environmental pollutant cadmium exposure, and identified KWMTBOMO12902 as a strong candidate gene for breeding applications in sericulture industry. Collectively, our results provide a novel and versatile approach for functional B. mori genomics, as well as a powerful resource for identifying the potential of key candidate genes for improving various economic traits. This study also shows the effectiveness, practicality, and convenience of large-scale mutant libraries in other nonmodel organisms.

Significance of the clock gene period in photoperiodism in larval development and production of diapause eggs in the silkworm Bombyx mori.

Many insects living in seasonal environments sense seasonal changes from photoperiod and appropriately regulate their development and physiological activities. Genetic researches have indicated the importance of a circadian clock system in photoperiodic time-measurement for photoperiodic regulations. However, most previous studies have focused on the effects on a single photoperiodic phenotype, without elucidating whether the circadian clock is involved in the core photoperiodic mechanism or only in the production of one target phenotype, such as diapause. Here, we focused on two different phenotypes in a bivoltine Kosetsu strain of the silkworm Bombyx mori, namely, embryonic diapause and larval development, and examined their photoperiodic responses and relationship to the circadian clock gene period. Photoperiod during the larval stage clearly influenced the induction of embryonic diapause and duration of larval development in the Kosetsu strain; short-day exposure leaded to the production of diapause eggs and shortened the larval duration. Genetic knockout of period inhibited the short-day-induced embryonic diapause. Conversely, in the period-knockout silkworms, the larval duration was shortened, but the photoperiodic difference was maintained. In conclusion, our results indicate that the period gene is not causally involved in the photoperiodic response of larval development, while that is essential for the short-day-induced embryonic diapause.

JH degradation pathway participates in hormonal regulation of larval development of Bombyx mori following λ-cyhalothrin exposure.

λ-Cyhalothrin (λ-cyh), a widely utilized pyrethroid insecticide, poses serious threats to non-target organisms due to its persistence nature in the environment. Exposure to low concentrations of λ-cyh has been observed to result in prolonged larval development in Bombyx mori, leading to substantial financial losses in sericulture. The present study was undertaken to elucidate the underlying mechanisms for prolonged development caused by λ-cyh (LC10) exposure. The results showed that the JH Ⅲ titer was significantly increased at 24 h of λ-cyh exposure, and the JH interacting genes Methoprene-tolerant 2, Steroid Receptor Co-activator, Krüppel-homolog 1, and JH binding proteins were also up-regulated. Although the target of rapamycin (Tor) genes were induced by λ-cyh, the biosynthesis of JH in the corpora allata was not promoted. Notably, 13 JH degradation genes were found to be significantly down-regulated in the midgut of B. mori. The mRNA levels and enzyme activity assays indicated that λ-cyh had inhibitory effects on JH esterase, JH epoxide hydrolase, and JH diol kinase (JHDK). Furthermore, the suppression of JHDK (KWMTBOMO01580) was further confirmed by both western blot and immunohistochemistry. This study has offered a comprehensive perspective on the mechanisms underlying the prolonged development caused by insecticides, and our results also hold significant implications for the safe production of sericulture.

SilkMeta: a comprehensive platform for sharing and exploiting pan-genomic and multi-omic silkworm data.

The silkworm Bombyx mori is a domesticated insect that serves as an animal model for research and agriculture. The silkworm super-pan-genome dataset, which we published last year, is a unique resource for the study of global genomic diversity and phenotype-genotype association. Here we present SilkMeta (http://silkmeta.org.cn), a comprehensive database covering the available silkworm pan-genome and multi-omics data. The database contains 1082 short-read genomes, 546 long-read assembled genomes, 1168 transcriptomes, 294 phenotype characterizations (phenome), tens of millions of variations (variome), 7253 long non-coding RNAs (lncRNAs), 18 717 full length transcripts and a set of population statistics. We have compiled publications on functional genomics research and genetic stock deciphering (mutant map). A range of bioinformatics tools is also provided for data visualization and retrieval. The large batch of omics data and tools were integrated in twelve functional modules that provide useful strategies and data for comparative and functional genomics research. The interactive bioinformatics platform SilkMeta will benefit not only the silkworm but also the insect biology communities.

Functional annotation of insecta transcriptomes: A cautionary tale from Lepidoptera.

Functional annotation is a critical step in the analysis of genomic data, as it provides insight into the function of individual genes and the pathways in which they participate. Currently, there is no consensus on the best computational approach for assigning functional annotation. This study compares three functional annotation methods (BLAST, eggNOG-Mapper, and InterProScan) in their ability to assign Gene Ontology terms in two species of Insecta with differing levels of annotation, Bombyx mori and Manduca sexta. The methods were compared for their annotation coverage, number of term assignments, term agreement and non-overlapping terms. Here we show that there are large discrepancies in gene ontology term assignment among the three computational methods, which could lead to confounding interpretations of data and non-comparable results. This study provide insight into the strengths and weaknesses of each computational method and highlight the need for more standardized methods of functional annotation.

Modeling the 3-dimensional structure of the silkworm's spinning apparatus in silk production.

The silk-spinning process of the silkworms transforms the liquid silk solution to a solid state under mild conditions, making it an attractive model for bioinspiration However, the precise mechanism behind silk expulsion remains largely unknown. Here we selected the silkworms as representative models to investigate the silk-spinning mechanism. We used serial block-face scanning electron microscopy (SBF-SEM) to reconstruct the three-dimensional structures of the spinnerets in silkworms at various stages and with different gene backgrounds. By comparing the musculature and duct deformation of these spinneret models during the spinning process, we were able to simulate the morphological changes of the spinneret. Based on the results, we proposed three essential factors for silkworm spinning: the pressure generated by the silk gland, the opening duct, and the pulling force generated by head movement. Understanding the silkworm spinning process provides insights into clarify the fluid-ejecting mechanism of a group of animals. Moreover, these findings are helpful to the development of biomimetic spinning device that mimics the push-and-pull dual-force system in silkworms. STATEMENT OF SIGNIFICANCE: The silkworms' spinning system produces fibers under mild conditions, making it an ideal candidate for bioinspiration. However, the mechanism of silk expulsion is unknown, and the three-dimensional structure of the spinneret is still uncertain. In this study, we reconstructed a detailed 3-dimensional model of the spinneret at near-nanometer resolution, and for the first time, we observed the changes that occur before and during the silk-spinning process. Our reconstructed models suggested that silkworms have the ability to control the spinning process by opening or closing the spinning duct. During the continuously spinning period, both the pressure generated by the silk gland and the pulling force resulting from head movement work in tandem to expel the silk solution. We believe that gaining a full understanding of the spinning process steps can advance our ability to spin synthetic fibers with properties comparable to those of native fibers by mimicking the natural spinning process.

Rapid production of chimeric silkworm/spider silk with improved mechanical properties by infection of nonpermissive Bombyx mori with recombinant AcMNPV harboring native-size of spidroin genes.

Spider silks with excellent mechanical properties attract more attention from scientists worldwide, and the dragline silk that serves as the framework of the spider's web is considered one of the strongest fibers. However, it is unfeasible for large-scale production of spider silk due to its highly territorial, cannibalistic, predatory, and solitary behavior. Herein, to alleviate some of these problems and explore aneasy way to produce spider fibers, we constructed recombinant baculovirus Autographa californica multiple nucleopolyhedrovirus (AcMNPV) simultaneously expressing Trichonephila clavipes native ampullate spidroin 2 (MaSp-G) and spidroin 1 (MaSp-C) driven by the promoters of silkworm fibroin genes, to infect the nonpermissive Bombyx mori larvae at the fifth instar. MaSp-G and MaSp-C were co-expressed in the posterior silk glands (PSGs) of infected silkworms and successfully secreted into the lumen of the silk gland for fibroin globule assembly. The integration of MaSp-G and MaSp-C into silkworm silk fibers significantly improved the mechanical properties of these chimeric silk fibers, especially the strength and extensibility, which may be caused by the increment of ß-sheet in the chimeric silkworm/spider silk fiber. These results demonstrated that silkworms could be developed as the nonpermissive heterologous host for the mass production of chimeric silkworm/spider silk fibers via the recombinant baculovirus AcMNPV.

Single-nucleus sequencing of silkworm larval midgut reveals the immune escape strategy of BmNPV in the midgut during the late stage of infection.

The midgut is an important barrier against microorganism invasion and proliferation, yet is the first tissue encountered when a baculovirus naturally invades the host. However, only limited knowledge is available how different midgut cell types contribute to the immune response and the clearance or promotion of viral infection. Here, single-nucleus RNA sequencing (snRNA seq) was employed to analyze the responses of various cell subpopulations in the silkworm larval midgut to B. mori nucleopolyhedrovirus (BmNPV) infection. We identified 22 distinct clusters representing enteroendocrine cells (EEs), enterocytes (ECs), intestinal stem cells (ISCs), Goblet cell-like and muscle cell types in the BmNPV-infected and uninfected silkworm larvae midgut at 72 h post infection. Further, our results revealed that the strategies for immune escape of BmNPV in the midgut at the late stage of infection include (1) inhibiting the response of antiviral pathways; (2) inhibiting the expression of antiviral host factors; (3) stimulating expression levels of genes promoting BmNPV replication. These findings suggest that the midgut, as the first line of defense against the invasion of the baculovirus, has dual characteristics of "resistance" and "tolerance". Our single-cell dataset reveals the diversity of silkworm larval midgut cells, and the transcriptome analysis provides insights into the interaction between host and virus infection at the single-cell level.